Universality in the diffusion of knots 
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We have evaluated a universal ratio between diffusion constants of the ring polymer with a given knot K and a 
linear polymer with the same molecular weight in solution through the Brownian dynamics under hydrodynamic 
interaction. The ratio is found to be constant with respect to the number of monomers, N, and hence the estimate 
at some N should be valid practically over a wide range of N for various polymer models. Interestingly, the 
ratio is determined by the average crossing number (Nac) of an ideal conformation of knotted curve K, i.e. 
that of the ideal knot. The Nac of ideal knots should therefore be fundamental in the dynamics of knots. 

PACS numbers: 82.35.Lr,05.40.Fb,05.20.-y 



I. INTRODUCTION 

Novel knotted structures of polymers have recently been 
found in various research fields such as DNA, proteins and 
synthetic polymers. 1,2,3,4 The topology of a ring polymer is 
conserved under thermal fluctuations in solution and repre- 
sented by a knot . 5 i 6 i 7 i 8 i 9 Topological constraints may lead to 
nontrivial statistical mechanical and dynamical properties of 
ring polymcn 6 - 10 - 11 - 12 - 13 - 14 - 15 - 16 - 17 - 18 - 19 . 

Recent progress in experiments of ring polymers should 
be quite remarkable. Diffusion constants of linear, relaxed 
circular and supercoiled DNAs have been measured quite 
accurately. 20 Here the DNA double lelices are unknotted. Fur- 
thermore, hydrodynamic radius of circular DNA has also been 
measured. 21 Ring polymers of large molecular weights are 
synthesized not only quite effectively 22 but also with small 
dispersions and of high purity j 2 ^ 4 - Circular DNAs with vari- 
ous knot types are derived, and they are separated into knot- 
ted species by gel electrophoresis. 25 We should remark that 
synthetic ring polymers with nontrivial knots have not been 
synthesized and separated experimentally, yet. However, it is 
highly expected that ring polymers of nontrivial knot types 
should be synthesized and their topological effects will be 
confirmed experimentally in near future. 

In the paper we discuss diffusion constant Dk of a ring 
polymer with fixed topology K in good solution for various 
knot types. We evaluate it numerically via the Brownian dy- 
namics with hydrodynamic interaction in which bond cross- 
ing is effectively prohibited through the finite extensible non- 
linear elongational (FENE) potential. 26 We evaluate diffusion 
constant Dl of a linear polymer with the same molecular 
weight, and derive ratio Dk/Dl. The ratio should corre- 
sponds to a universal amplitude ratio of critical phenomena 
and play a significant role in the dynamics of knotted ring 
polymers. According to the renormalization group arguments, 
ratio Dk/Dl should be universal if the number of monomers, 
N, is large enough. 27 i 28 i 29 

The ratio Dk/Dl may have some experimental applica- 
tions. Ring polymers of different knot types can be separated 
experimentally with respect to their topologies by making use 
of the difference among the sedimentation coefficients, which 
can be calculated from the diffusion constants.— Here we re- 
mark that the diffusion constant of a ring polymer under no 



topological constraint , Dr, and that of the corresponding 
linear polymer has been numerically evaluated, and the ratio 
C = D r /Dl has been studied.^ - 3 ^ 

Through simulation we find that ratio Dk/Dl is approxi- 
mately constant with respect to N for various knots. Thus, if 
we evaluate ratio Dk/Dl at some value of TV, it is practically 
valid for other values of N. We can therefore predict the dif- 
fusion constant Dk of a polymer model at some value of N, 
multiplying the ratio Dk/Dl by the estimate of Dl- Here 
we remark that the value of Dl may depend on the number N 
and on some details of the model j^ 3 -* 3 - 4 - 

Furthermore, we show numerically that ratio Dk/Dl is 
a linear function of the average crossing number (Nac) of 
the ideal knot of K, an ideal configuration of knotted curve 
K, which will be defined shortly. Since the ratio Dk/Dl 
is almost independent of N, it follows that the linear fitting 
formula should be valid practically in a wide range of finite 
values of N. Thus, the ideal knot of a knotted curve K should 
play a fundamental role in the dynamics of finite- size knotted 
ring polymers in solution. 

Let us introduce the ideal knot, briefly. For a given knot K 
it is given by the trajectory that allows maximal radial expan- 
sion of a virtual tube of uniform diameter centered around the 
axial trajectory of the knot K&& We define the Nac of a 
knotted curve as follows: We take its projection onto a plane, 
and enumerate the number of crossings in the knot diagram 
on the plane. Then, we consider a large number of projections 
onto planes whose normal vectors are uniformly distributed 
on the sphere of unit radius, and take the average of the cross- 
ing number (Nac) over all the normal directions. 

The paper consists of the followig: In section II, the sim- 
ulation method is explained. In section III, we present the 
estimates of the diffusion constant of a ring polymer in solu- 
tion of knot type K for various knot types. Then, we show 
numerically that the graph of Dk/Dq is almost independent 
of N, and also that ratio Dk/Dl is fitted by a linear function 
of Nac of the ideal knot of K. We also discuss the simula- 
tion result in terms of the ratio of equivalent radii, 37 clg/clt, 
which corresponds to the universal ratio of the radius of gyra- 
tion to the hydrodynamic radius. 33 We shall define the equiva- 
lent radii explicitly in section III. Finally, we give conclusion 
in section IV. 

Throughout the paper, we employ the symbols of knots fol- 



lowing Rolfsen's textbook 38 , as shown in Figure 1. 



II. SIMULATION METHOD 

The ring polymer is modeled as a cyclic bead- and- spring 
chain with N beads connected by N FENE springs with force 
given by 



knot type 



D 



(R 2 g) 



F(r) = -Hr/(l-r 2 /r 2 max ) 



(1) 



where r = \r\. We denote by b the unit of distance, which 
gives the average distance between neighboring monomers 
approximately. We set constants H and r max by H = 
30/c#T/6 2 and r ma x = 1-36. We assume the Lennard- Jones 
(LJ) potential by 



V(rij) = 4e L j{(cr LJ /r ij ) 12 - (cr L j/r^) } . 



(2) 



Here r^ is the distance of beads i and j, and 6lj and ctlj 
denote the minimum energy and the zero energy distance, 
respectively. 39 We set the Lennard- Jones parameters as gl j = 
0.8b and clj = 0.1 k&T so that they give good solvent 
conditions. 40 Here ks denotes the Boltzmann constant. 

We employ the predictor-corrector version 41 of the Ermak- 
McCammon algorithm for generating time evolution of a 
ring polymer in solution. The hydrodynamic interaction 
is taken into account through the Rotne-Prager-Yamakawa 
tenso r 39 i 42 i 43 where the bead friction is given by ( = 67rr] s a 
with the bead radius a = 0.2576 and a dimensionless hydro- 
dynamic interaction parameter h* = (£/6irri s )<\/H/irkBT = 
0.25. 

In the present simulation, physical quantities are given in 
dimensionless units such as in Ref. 39 . We divide length by 
b, energy by ksT and time by (b 2 /k^T. Let us indicate di- 
mensionless quantities by an asterisk as superscript. We have 
H* = 30, r^ ax = 1.3. We take the simulation time step 
At* = 10- 4 . 

We have set the FENE potential so that the topology of 
the ring polymer should be effectively conserved. 26 However, 
bond crossing may occur with a very small probability. Cal- 
culating knot invariants, we have confirmed that the fraction 
of nontrivial knots is very small. If the initial knot type is the 
trivial knot it is given by 10 -8 ~ 10 -7 , and if the initial knot 
type is a nontrivial knot, it is given by approximately 10 -7 . 



III. SIMULATION RESULTS 



We define the diffusion constant of a polymer by 



D= liml((r G (i)-r G (0)) 2 ). 



(3) 



Here rc(t) denotes the position vector of the center of mass 
of the polymer at time t. Making use of 0]) we have evaluated 
diffusion constants Dl and Dk- 

The estimates of diffusion constants Dl and Dk at N = 45 
are listed in Table 1 together with those of the mean square 



linear 


0.12038 ±0.00085 9.33029 ± 0.03219 





0.13059 ±0.00089 5.26539 ± 0.00913 


3i 


0.14530 ±0.00079 3.21052 ± 0.00505 


4i 


0.14876 ±0.00074 2.78817 ± 0.00160 


5i 


0.15277 ±0.00085 2.72300 ± 0.00261 


5 2 


0.15640 ±0.00078 2.61427 ± 0.00132 


6i 


0.15927 ±0.00078 2.47449 ± 0.00137 


6 2 


0.15902 ±0.00095 2.37272 ± 0.00126 


7i 


0.16416 ±0.00083 2.47162 ± 0.00109 



TABLE I: Estimates of diffusion constants Dl and Dk, and the 
mean square radius of gyration (R%) and for a linear polymer of 
N — 45 and ring polymers of N — 45 with various knot types. 



radius of gyration (Rq). The data of Dl and Dk are plotted 
against N in Figure 2. The fitting curves to them are given 
by D = aN~ v (l ± bN~ A ). Here the errors of the diffusion 
constants are as small as 10 -4 . 

Ratios Dk/Dl should correspond to universal amplitude 
ratios in critical phenomena. Numerically we find that ratio 
Dk ± /Dk 2 of two different knots K\ and K 2 is almost con- 
stant with respect to N, at least in the range investigated . For 
instance, the graph of ratio Ds ± /Do versus N and that of ra- 
tio D^ ± /Do versus N for the data are almost flat, as shown 
in Figures 3 and 4, respectively. Here 0, 3i and 4i denote 
the trivial, the trefoil and the figure-eight knot, respectively, 
as shown in Figure 1. The numerical values of D^/Do are 
given from 1.14 to 1.17 in Figure 3, and those of D^ ± /Do are 
given from 1.14 to 1.21 in Figure 4. Thus, the estimate of 
Dk/Do evaluated at a value of N, say N = 45, for some 
knot K should also be valid at other finite values of N, since 
it is almost independent of N. 

For the diffusion constant of a ring polymer, Dr, the ra- 
tio Dr/Dl should correspond to a universal amplitude ratio 
and should be universal if N is large enough i 27 i 28 i 29 For the 
diffusion constant Dr, there is no topological constraint in 
the ring polymer model and Dr does not mean Dk of a knot 
if i 26 i 3Q i 31 In the previous simulation 26 it has been found that 
ratio Do/Dl is given by about 1.1 for the present polymer 
model and almost independent of N within the range investi- 
gated. 

From the numerical observations and the RG arguments, we 
have two conjectures: (A) Do/Dl should be given by 1.1 for 
some wide range of finite values of N and also in the large N 
limit; (B) ratio Dk/Do for a nontrivial knot K should remain 
almost the same value in a wide range of finite values of N, 
i.e. the TV-dependence should be very small. 

Quite interestingly we find that ratio Dk/Dl can be ap- 
proximated by a linear function of the average crossing num- 
ber (Nac) of ideal knots, i.e. the ideal representations of the 
corresponding knots. In Figure 5 simulation data of Dk/Dl 
are plotted against Nac of ideal knots. We find that the data 
points are fitted well by the following empirical formula: 



Dk/Dl = a + bN AC . 



(4) 



Here, the estimates of a and b are given in the caption of Fig- 
ure 5. Thus, the diffusion constant Dk of a knot K can be 
estimated in terms of the Nac of the ideal knot of K. 

Let us discuss the \ 2 values. We have \ 2 — 2 for the fitting 
line of Figure 5, which is for the data of N = 45. For the 
data of N = 36 we have a good fitting line with x 2 = 3. 
The estimates of a and b for N = 36 are similar to those for 
N = 45. Thus, we may conclude that the graph of Dk/Dl 
versus Nac is fitted by a linear line. 

For a finite value of N, we can estimate the diffusion con- 
stant Dk of a knot K through formula (0]) by the Nac of the 
ideal knot of K. Here we have assumed that coefficients a 
and b of © are independent of N since the graphs of Dk/Dq 
and Dq/Dl are almost flat with respect to N. In fact, there 
is almost no numerical support for suggesting a possible TV- 
dependence of a and b, directly. 

We thus summarize the simulation results so far as follows: 
ratio Dk/Dq for a knot K should be almost constant with re- 
spect to N in a wide range of N and can be expressed by the 
linear function of Nac of ideal knots. Eq. Q should be use- 
ful in separating synthetic ring polymers into various knotted 
species by making use of the difference among sedimentation 
coefficients. 

Ideal knots should play a fundamental role in the dynam- 
ics of knotted ring polymers in solution. In fact, we have 
shown it for the diffusion constants. In experiments of gel 
electrophoresis drift velocities of different knots formed on 
the same DNA molecules were shown to be simply related 
to the Nac of ideal knots. 25 The two independent results sug- 
gest the importance of the Nac of ideal knots in the dynamics 
of knotted ring polymers, although the physical situations are 
different. 

Let us now discuss the simulation results from the view- 
point of equivalent radii. 37 The equivalent radius for any so- 
lution property is the radius of a spherical particle having the 
same value of solution property as that of the macromolecule 
under consideration. The ratio of equivalent radii should be 
universal, and it should play a similar role as the universal 
amplitude ratio such as the ratio of diffusion constants. 33 We 
define equivalent radii ac and ar explicitly by 



a G 



\{Rl), 



k B T 

67TT] S D 



(5) 
(6) 



Here ac and cll corresponds to the radius of gyration Rq = 
-\/(Rg) and the translational friction coefficient D, resepc- 
tiveley. The ratio ac/ar corresponds to the ratio of the radius 
of gyration to the hydrodynamic radius, and should be univer- 
sal. 

The numerical estimates of ac/ar for N = 45 for the 
present simulation are listed in Table 2 for linear and ring 
polymers with various knot types. In Figure 6, the ratio 
clg/clt is plotted against the number of segments, N, for lin- 



ear and ring polymers with various knot types. Interestingly, 

the graphs show a weak 7V-depndence. They are fitted by a 

function ac/ar = a(l — bN~ c ) with parameters a, b and c 

knot type clg/clt 

linear 1.8471 ± 0.00475 

1.5053 ±0.01869 

3i 1.3079 ±0.01151 

4i 1.20824 ±0.00055 

5i 1.26636 ±0.00077 

5 2 1.27029 ±0.00052 

6i 1.25865 ±0.00084 

6 2 1.23046 ±0.00058 

7i 1.29643 ±0.00050 



TABLE II: Estimates of ratio a G /a T = y/b(R%)/3D*/a* for a 
linear polymer of N = 45 and ring polymers of N — 45 with various 
knot types. Here, (R%) = (R G )/b 2 , D* = 67Tf] s aD/k B T and 

a* = a/b. 



being positive. It suggests that the graphs become constant 
with respect to N if N is large enough. We thus expect that 
the ratio clg/clt in the large TV limit should be universal. 

It is interesting to note in Figure 6 that the estimate of ratio 
ac/dr in the large TV limit is distinct for the different topolo- 
gies such as linear polymers and ring polymers of the trivial 
and trefoil knots. The ratio could thus be useful for detecting 
the knot type of a ring polymer in solution. 



IV. CONCLUSION 

We have evaluated universal ratios among the diffusion 
constants of knotted ring polymers in good solution for sev- 
eral knots, where bond crossing is effectively prohibited in 
the Brownian dynamics under hydrodynamic interaction. The 
universal ratio of diffusion constants Dk/Dl is almost con- 
stant with respect to the number of polymer segments, N. 
Moreover, it is found that the ratio Dk/Dl is determined by 
the Nac of the ideal knot of K. Through the linear relation, 
we can estimate the diffusion constant of a given knot. 
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FIG. 1 : Figures of a linear polymer and knotted ring polymers with 
the symbols of knots given in Rolfsen's textbook. 
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FIG. 2: Diffusion constants of linear and knotted ring chains with 
knots 0, 3i, 4i, 5i, 6i and 7i, versus N. Fitted by D = aN~ u (l + 
bN~ A ) with the following best estimates: For a linear chain, a — 
0.90 ± 0.23, v = 0.53 ± 0.06, b = 0.51 ± 0.93, A = 1.14 ± 2.39, 
X 2 = 17; for the trivial knot (0), a = 1.03 ± 1.11, z/ = 0.55 ± 0.18, 
6 = 0.14±0.78, A = 0.60=b6.09, \ 2 = 28; for the trefoil knot (3 1) 
a = 1.00±3.87, v = 0.52±0.67, b = 1.18±1.12, A = 0.77±6.09, 
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FIG. 3: Ratio Dz x /D of diffusion constants for the trefoil knot (3i) 
and the trivial knot (0) versus the number of segments N. Fitting 
curve is given by D 3 JD = a(l + 6iV" c ), where a = 1.07 ±0.64, 
b = 0.25 ± 0.39, and c = 0.39 ± 3.29 with x 2 = 6. 
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FIG. 4: Ratio D^ /Do of diffusion constants for the figure-eight 
knot (4i) and the trivial knot (0) versus the number of segments N. 
Fitting curve is given by D^ /D = a(l-\-bN~ c ) where a = 1.02± 
0.56, b = 1.76 ± 8.26, and c = 0.70 ± 2.58 with x 2 = 0.03. 
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FIG. 5: Dk/Dl versus the average crossing number (Nac) of ideal 
knot K for N = 45: The data are approximated by Dk/Dl = 
a + bN A c where a = 1.11 ± 0.02 and b = 0.0215 ± 0.0003 with 
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FIG. 6: clg/clt of linear and knotted ring chains with knots 0, 3i, 
4i, 5i, 6i and 7i, versus AT. Fitted by aG/ar = a(l — bN~ c ) with 
the following best estimates: For a linear chain, a = 2.37 ± 0.29, 
6 = 0.52 ± 0.03, c = 0.22 ± 0.10, x 2 = 55; for the trivial knot (0), 
a = 1.51±0.03, b = 0.72±0.23, c = 0.89±0.22, X 2 = 14; for the 
trefoil knot (3i), a = 1.30±0.03, b = 2.10±3.71, c = 1.27±0.76, 
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